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ABSTRACT 

Context. Massive galaxy clusters at intermediate redshift can magnify the flux of distant background sources by several magnitudes. 
Aims. We exploit this effect to search for lensed distant supernovae that may otherwise be too faint to be detected. 
Methods. A supernova search was conducted at near infrared wavelengths using the ISAAC instrument at the VLT. The massive galaxy 
clusters Abell 1689, Abell 1835, and AC1 14 were observed for a total of 20 hours to search for supernovae in gravitationally magnified 
background galaxies. The observations were split into individual epochs of 2 hours of exposure time, separated by approximately one 
month. Image-subtraction techniques were used to search for transient objects with light curve properties consistent with supernovae, 
both in our new and archival ISAAC/VLT data. The limiting magnitude of the individual epochs was estimated by adding artificial 
stars to the subtracted images. Most of the epochs reach 90% detection efficiency at SZ( J) = 23.8 - 24.0 mag (Vega). 
Results. Two transient objects, both in archival images of Abell 1689 and AC1 14, were detected. The transient in AC1 14 coincides - 
within the position uncertainty - with an X-ray source and is likely to be a variable AGN at the cluster redshift. The transient in Abell 
1689 was found at SZ = 23.24 mag, ~ 0.5"away from a galaxy with photometric redshift z ga i = 0.6 + 0.15. The light curves and the 
colors of the transient are consistent with a reddened Type IIP supernova at redshift z = 0.59 ± 0.05. The lensing model of Abell 1689 
predicts ~1.4 mag of magnification at the position of the transient, making it the most magnified supernova ever found and only the 
second supernova found behind a galaxy cluster. 

Conclusions. Our pilot survey has demonstrated the feasibility to find distant gravitationally magnified supernovae behind massive 
galaxy clusters. One likely supernova was found behind Abell 1689, in accordance with the expectations for this survey, as shown in 
an accompanying analysis paper. 
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1. Introduction 

Gravitational lensing by massive galaxy clusters provides unique 
opportunities to study extremely distant background objects. 
The modeling of the most massive galaxy clusters at interme- 
diate redshifts shows that the flux of background sources with 
lines-of-sight close to the critical lines can be magnified by 
more than 3 magnitudes. These powerful natural gravitational 
telescopes have been successfully used to identify and study 
gravitationally-lensed galaxies at very high redshift (e.g., Franx 
et al., 1997; Ellis et al., 2001; Hu et al., 2002; Kneib et al., 2004; 
Richard et al., 2006, 2008). 
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Repeated observations of massive galaxy clusters may also 
allow discovery of very distant supernovae (SNe) at redshifts 
beyond z - 1.5-2, which may be too faint to be detected with- 
out the lensing magnification (e.g., Kovner & Paczynski, 1988; 
Kolatt & Bartelmann, 1998; Saini et al., 2000; Sullivan et al., 
2000; Gal- Yam et al., 2002; Gunnarsson & Goobar, 2003). Gal- 
Yam et al. (2002) have searched repeated observations of galaxy 
clusters from the HS T archive and found several SNe, but only 
one of them in a background galaxy, which has a redshift close 
to that of the lensing cluster. 

Detecting and studying SNe at very high redshifts is impor- 
tant for several reasons: 

1. The progenitors of core-collapse (CC) SNe are massive 
short-lived stars. The CC SN rate thus reflects the ongoing star 
formation rate (SFR) and can be used to measure it. The first 
results on the SFR derived from CC SNe (Dahlen et al., 2004) 
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show an increase in the range z ~ 0.3 - 0.7, consistent with the 
measurements from galaxy luminosity densities. The magnifi- 
cation provided by galaxy clusters may help to detect a popu- 
lation of CC SNe at z > 1 for the first time and to obtain an 
independent measurement of the cosmic SFR in a redshift range 
where current theoretical predictions do not give consistent re- 
sults (e.g., Dahlen et al., 2004; Mannucci et al., 2007; Kobayashi 
& Nomoto, 2008). It may also be possible to put constraints on 
the rates of the very luminous Type Iln and pair-production SNe, 
and on Type Ic hypernovae. In addition, there is growing evi- 
dence that Type IIP SNe could be calibrated as standard candles 
(Hamuy & Pinto, 2002; Nugent et al., 2006; Poznanski et al., 
2009; Olivares, 2008) and may turn out to be more useful cos- 
mological probes than Type la SNe (SNe la) at high redshift 
because SNe la may be very rare or even absent at z > 2-3 
(Kobayashi et al., 1998). 

2. Type la SNe have been used to measure the cosmological 
parameters; however, the scenario(s) leading to the explosion of 
SNe la is still unknown. Models predict that different progenitor 
scenarios should have different delay times r between the on- 
set of star formation and the SN la explosion. By measuring the 
SN la rate and comparing with the SFR, it is possible to set con- 
straints on r, hence on progenitor scenarios (see e.g., Strolger 
et al., 2004). The predicted SN la rate is most sensitive to the 
delay time for z > 1.5, but, as for CC SNe, different SN la rate 
predictions significantly diverge at z > 1.5 (Kobayashi et al., 
1998; Dahlen et al., 2008; Mannucci et al., 2007; Kobayashi 
& Nomoto, 2008; Neill et al., 2006; Scannapieco & Bildsten, 
2005). Strolger et al. (2004) and Dahlen et al. (2008) find a long 
delay time t ~ 3 - 4 Gyr from the low number of SNe found at 
z > 1 .4. This implies that SN la progenitors are old, low-mass 
stars and that there should be a steep decline in the SN la rate 
at z > 1.5 (as in the model of Kobayashi et al. 1998). However, 
Mannucci et al. (2005) and Sullivan et al. (2006) present evi- 
dence for a dependence of the SN la rate on the SFR, thus re- 
lating a significant fraction of SNe la to young massive stellar 
populations. It is expected then that the SN la rate will increase 
at z > 1.5, although Mannucci et al. (2007) suggest that a con- 
siderable fraction of SNe la may be missed because of dust ex- 
tinction. Measurement of the SN la rate at z > 1.5 will thus be 
critical for distinguishing between the models. 

3. Modeling of the mass distribution of the best-studied clus- 
ters shows that the magnification of lensed SNe can be estimated 
to an accuracy of < 0.25 mag or better depending on the SN lo- 
cation and the accuracy of the mass model (Kneib et al., 1996; 
Saini et al., 2000). Thus, lensed SNe can be used as distance 
indicators and significantly increase the leverage of the current 
SN la Hubble diagram. 

4. The time delay between multiple images of highly magni- 
fied galaxies could be as shotr as weeks or months. Detection of 
multiple SN images is then possible, and the time delay between 
them could be used to constrain the Hubble parameter (Refsdal, 
1964). Because of the transient nature of SNe, the uncertainty in 
the time delay will be small (Goobar et al., 2002; Mortsell et al., 
2005). 

In this paper we describe a pilot near-infrared (NIR) survey 
designed to find gravitationally magnified SNe in galaxies be- 
hind the massive, well-studied galaxy clusters Abell 1689, Abell 
1835 and AC1 14. Here, we describe the observing strategy, data 
reduction, and the SN search efficiency of the survey. In an ac- 
companying paper (Goobar et al., 2009, hereafter Paper II) we 
present a detailed analysis of the survey and the implications for 
future deeper and wider surveys. A feasibility study of the po- 
tential for improving the mass models of clusters of galaxies and 



measuring cosmological parameters using lensed SNe will be 
presented in Riehm et al. (in preparation, Paper III). Throughout 
the paper we assume the concordance cosmological model with 
Qm = 0.3, Qa = 0.7, w = -1 and h = 0.7, and we use magni- 
tudes in the Vega system. 



2. Observations 

For this pilot survey, we selected three of the best-studied mas- 
sive galaxy clusters at intermediate redshift: Abell 1689 (z = 
0.183), Abell 1835 (z = 0.253), and AC114, also known as 
Abell S 1077 (z — 0.312). The mass distributions of these clusters 
have been extensively modeled (Broadhurst et al., 2005; Richard 
et al., 2006; Limousin et al., 2007) and the magnification maps 
can be computed as a function of source redshift. Deep multi- 
color optical and NIR observations are available for all clusters, 
including HS T data, making it possible to compute photometric 
redshifts of most of the galaxies in the FOV. In addition, spec- 
troscopic redshifts of multiple images of many strongly lensed 
background galaxies are available. 

At the targeted redshifts, beyond z ~ 1 , the U V/optical wave- 
lengths, where most of the SN luminosity is emitted are red- 
shifted to the near infrared (NIR) domain, and the search should 
be ideally done in the NIR. Even when including the lensing 
magnification, supernovae at z ~ 1 are expected to be fainter 
than 23 mag in the NIR, and thus 8-m class telescope observa- 
tions are needed. For the observations described in this paper, we 
used VLT/ISAAC (Moorwood et al., 1998). The Einstein radii 
of the three clusters at high redshift are ~ 30 - 50", and the 
large-magnification regions fit well into the ISAAC 2.5'x2.5' 
field-of-view (FOV). In addition, archival NIR observations at 
VLT/ISAAC under excellent seeing conditions are available for 
all clusters and can be used as reference images for the SN 
search. 

To test the feasibility of the project, we used archival 
VLT/ISAAC observations of Abell 1835 obtained with an SZ 
filter 1 (Richard et al., 2006). The observations were obtained in 
two runs one month apart and with total exposure times of ~ 2 
and ~ 4 hours. The two combined images were subtracted and 
searched for transient sources. None were found, but these obser- 
vations allowed us to test the data reduction and the image sub- 
traction scheme on real data, and to realistically estimate the SN 
survey depth. By adding artificial stars we could establish that 
we would have a good detection efficiency down to SZ ^ 24.5 
mag. Using SZ 24.5 mag as the projected photometric depth 
of the survey, the cluster magnification map from Richard et al. 
(2006) and the absolute magnitudes of the different SN types, a 
first estimate of the number of SNe that could be detected in our 
survey was obtained following the procedure in Paper II. The 
results were encouraging and we estimated that a few SNe per 
cluster per year could be detected. For the highest redshifts, be- 
sides Type la SNe, the survey should be most sensitive to Type 
Iln and IIL SNe, which have very blue spectral energy distri- 
bution (SED) around maximum light. Close to peak brightness, 
Type IIP SNe should also be within reach. Since significant star 
formation may take place in the blue arcs of the strongly lensed 
galaxies behind the clusters, several magnified Type II SNe are 
expected in our survey, possibly with multiple images within our 
observational window. 



1 The SZ filter is a broadband filter with central wavelength 1.06/jm 
and width 0.13yum. It is most often used as an order-sorting filter for 
ISAAC NIR spectroscopy but is also offered as an imaging filter. 
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After the preliminary study, we decided to observe our tar- 
gets with 2 hours of total exposure time per observation with 30 
days cadence over the period when the clusters were observable. 
The search was thus performed in a 'rolling search' scheme, 
where new SNe are discovered and old SNe are followed us- 
ing the same observations. Abell 1689 and Abell 1835 were ob- 
served with the SZ filter and AC 114 with J to be able to use 
the archival VLT/ISAAC observations as references. Details of 
our observations are given in Table 1 (the 2007 data), along with 
other NIR and optical data from the ESO archive. No spectral 
or multi-color observations were obtained and for this pilot run, 
the SN identification is done using the information on light curve 
shapes, absolute magnitudes and host galaxy type and redshift, 
if known. 

3. Data reduction 

Because of the presence of bright galaxies with extended ha- 
los at the cluster centers, the best observing strategy is beam- 
switching. However, this observing mode doubles the execution 
time needed for the observations. We therefore used a simple 
dithering scheme without beam-switching for separate sky mea- 
surements. The exposures were quasi-randomly dithered within 
a box of size 40"- larger than the typical size of the bright cluster 
galaxies - and we adopt a slightly revised version the reduction 
procedure described by Richard et al. (2006) in order to opti- 
mize the background subtraction and the detectability of faint 
point sources. 

All images were first corrected for the 'electrical ghost' 
effect with the recipe ghost within the ESO software library 
Eclipse 2 . The effect consists of an additional signal, which on 
a given row is proportional to the sum of the intensity along this 
row and the row 512 rows away. It has a reproducible behavior 
and can be removed with a simple algorithm 3 . Dark current and 
flat-field corrections were applied using observations obtained 
as a part of the ISAAC standard calibration plan. The most chal- 
lenging part of the NIR data reduction is the background subtrac- 
tion. We performed this with the XDIMSUM package in IRAF 4 
in a 3-step process: 

1. For each image, the 8 exposures closest in time were scaled 
to the same mode and median combined to estimate the sky 
background. The combined image was scaled to the mode of 
the image from which the background is to be removed and 
was subtracted from it. The sky-subtracted images were then 
combined with integer pixel shifts to preserve the character- 
istics of the noise. The ISAAC pedestal is highly variable in 
time and is very difficult to remove accurately. In some of 
the images, significant residuals were found along the bot- 
tom and the middle rows of the array, where the readout of 
the two halves of the array starts. It was then necessary to 
correct for this effect before the final combination. At the 
first step this was done by simply taking a median along the 
image rows with 2<x clipping and subtracting it. 

2. The combined image from Step 1 was used to detect ob- 
jects and create an object mask. For this XDIMSUM uses 

2 available at http : //www. eso . org/projects/aot/eclipse/ 

3 see http : //www. eso . org/sci/facilities/paranal/ 
instruments/ i saac/ doc/ drg/html/node3 3 . html 

4 All data reduction and calibration was done in IRAF and with 
our own programs written in IDL. IRAF is distributed by the 
National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 



Table 1. Data used to search for transient objects. 



Date Exposure Seeing 90% detection 

[min] [arcsec] efficiency [mag] 

Abell 1689 - VLT/ISAAC SZ-band 



pointing 1 



2003 02 09 


159 


0.52 24.28, transient 


2003 04 27 


43 


0.43 24.28, transient 


2004 01 13 


43 


0.52 23.58, non-detect 


2004 02 14 


43 


0.58 23.64, non-detect 




pointing 2 


2003 01 16 


43 


0.58 23.48 


2003 02 15 


43 


0.50 23.60 


2003 04 27 


86 


0.44 23.64 


2004 01 12 


43 


0.55 23.64 




pointin 


g to cluster core 


2007 04 08 


117 


0.64 23.95 


2007 05 14/15 


117 


0.65 23.95 


2007 06 06 


39 


0.70 23.15 


AC 1 14 - VLT/ISAAC 7-band 


2002 08 20 


108 


0.49 23.87 


2007 07 13* 


234 


0.43 24.04 


2007 08 09 


117 


0.73 23.79 


2007 09 02 


117 


0.55 23.83 


2007 09 28 


117 


0.46 24.04 


Abell 1835 - 


VLT/ISAAC 5 Z-band 


area 1 






2004 04 20 


231 


0.49 24.06 


2004 05 15 


135 


0.62 24.06 


2007 04 18 


117 


0.79 23.80 


2007 05 18 


78 


0.74 23.83 


2007 07 18 


117 


0.62 23.80 


area 2 






2007 04 18 


117 


0.79 23.70 


2007 05 14/18 


60 


0.80 23.45 


2007 07 18 


117 


0.62 23.70 



" - average of observations obtained on 5, 1 1 and 15 February; 
6 - average of observations obtained on 1 1,12,13 and 15 July. 



a simple method where all pixels above a given threshold, 
typically 2 times the standard deviation of the sky, are as- 
sumed to belong to objects. The mask is then deregistered 
to each individual image. The same operations as in Step 1 
were then repeated, but now all pixels belonging to objects 
were excluded when performing the background estimation 
and residual pedestal subtraction. The only difference was 
that, instead of the median, the mode along the image rows 
was used to remove the residual pedestal. 
3. The combined image from Step 2 was used to create an im- 
proved object mask using the more sophisticated OBJMASK 
task. This task detects connected pixels and can also apply 
object growing, which in our case was important for masking 
the extended halos of the cluster galaxies. The background 
subtraction was repeated with the improved object mask to 
arrive at our final sky-subtracted images. 

For the final combination, we assigned weights to the indi- 
vidual images in order to optimize the detection of faint point 
sources. The weights are inversely proportional to the product of 
the square of the seeing, the variance of the sky background and 
the multiplicative factor that brings the images to the same flux 
scale. The archival ISAAC observations, which we used as refer- 
ence images and when possible to search for SNe, were reduced 
following the same recipe. 

The photometric calibration of the images was obtained 
with standard stars from Persson et al. (1998) and a few from 
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Magnitude 

Fig. 2. The subtraction of the 7-band images of AO 14. The arrow 



Fig. 1. The detection efficiency as a function of magnitude for the ob- shows the transient. The ACS z-band image is also shown, 
servations of AC 114. The results for Abell 1835 and Abell 1689 are 
qualitatively the same and are not shown. 



Hawarden et al. (2001) observed in each photometric night as 
part of the ISAAC standard calibration plan. Because the stan- 
dard stars do not have calibrated SZ magnitudes, we used J- 
band magnitudes to calibrate the SZ observations of Abell 1835 
and Abell 1689. The majority of the Persson et al. (1998) stars 
have J - H ss 0.3 mag, which corresponds to late F to early G 
type stars. Synthetic SZ and / photometry was computed using 
spetrophotometry of F8V to G5V stars from Pickles (1998), a 
synthetic spectrum of the Sun from Kurucz web site 5 and the 
Vega spectrum from Bohlin (2007), assuming all colors of Vega 
to be zero. The computed synthetic SZ - J colors ranged be- 
tween 0.1 and 0.2 mag, and we corrected the zero points esti- 
mated from the / magnitudes of the Persson et al. (1998) stars 
by +0.15 mag to derive the SZ zero points. A correction of the 
same order was also estimated based on synthetic photometry of 
elliptical galaxies, as described in Richard et al. (2006). 

4. The transient search and survey depth 

To search for transient objects, the reference images were sub- 
tracted from our new images. The subtraction was done with 
Alard's optimal image subtraction software (Alard & Lupton, 
1998; Alard, 2000), slightly modified and kindly made avail- 
able to us by B. Schmidt. The reference image is geometrically 
aligned to the new image, convolved by a kernel such that the 
two images have the same PSF, and then scaled and subtracted. 
The residual images were searched for transient sources with 
Sextractor v2.5 (Bertin & Arnouts, 1996). The images were first 
convolved with a Gaussian filter with FWHM equal to the see- 
ing and then searched for sources consisting of at least four con- 
nected pixels lcr above the background. 

Because we expect most SNe to be very faint and close to 
the detection limit of the survey, we required the object to be 
detected in at least two epochs for a secure detection. Since SNe 
could have been recorded on the archival data that were used 
as reference, we also searched for negative peaks. Part of the 
archival observations were obtained during many separate runs 
and were also searched for transients. 

The detected sources were carefully examined visually and 
obvious artifacts and clear spurious detections near the underex- 
posed image edges discarded. The cores of the compact, bright 

5 available at http : / /kurucz . harvard . edu 



cluster galaxies could often not be subtracted cleanly, so sig- 
nificant residuals were left. Therefore, all "detections" close to 
badly subtracted cluster galaxies were rejected. After these cuts, 
only a few detections per image were typically left for closer 
examination. None of these potential transients was detected in 
more than one epoch. The careful visual examination (includ- 
ing blinking of the subtracted and the original images) showed 
that all but two candidates could not be visually distinguished 
from the background noise. That none of the other potential can- 
didates were picked up by Sextractor in another epoch makes 
us conclude that there was no evidence that any of these were 
real transient detections. The remaining two detections will be 
discussed in detail in the next section. 

To estimate the detection limits of the survey, for each epoch 
the PSF was computed using stars in the fields and 50 artificial 
stars with magnitudes between 22 and 25 were randomly added 
to the subtracted images. For each epoch, this was repeated 200 
times (10000 artificial stars), and the artificial stars were recov- 
ered with the same Sextractor parameters used to search for tran- 
sients. The detection efficiency was then computed in 0.2 mag 
bins as the fraction between the number of recovered stars and 
the number of added stars. Because we expect to find SNe on 
faint galaxies (if visible at all) with surface brightnesses close 
to the background, we only considered the artificial stars away 
from bright objects in the field. In addition, because we expect 
to find rather faint SNe close to the detection limit, we only con- 
sider the central parts of the images that received at least 90% of 
the total exposure time. In Table 1, we show the 90% detection 
limits and in Fig. 1, the detection efficiency curves for AC114 
are shown as an example. 

We note several complications. For AC 114 our first epoch 
had better seeing and photometric depth than the archival obser- 
vations. We therefore used our first epoch as the main reference 
image for the search. A Type IIP SN that had nearly constant 
brightness during our observational window could still be de- 
tected as a negative source in the subtraction of the archival ob- 
servations. 

The archival observations of Abell 1689 have been obtained 
at 4 pointings in order to cover the whole cluster. Each pointing 
covers ~ 1/4 of our new images and so each quarter was then 
treated separately using the respective reference image. The four 
combined reference images have nearly the same seeing, and the 
detection efficiencies were nearly the same. In Table 1 we give 
the average detection limits. 
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Table 2. Measured magnitudes of the transient in Abell 1689. 



Date 


Magnitude 


Filter 


Telescope/ 




[VegaJ 




Instrument 


2003 02 06 


24.09 ± 0.20 


/ 


VLT/FOSR2 


2003-02-9/10 


23.93 ± 0.08 


c 


VLT/FOSR2 


2003-02-26/27 


23.94 ± 0.09 


z 


VLT/FOSR2 


2003 02 09 


23.24 ± 0.08 


SZ 


VLT/ISAAC 


2003 04 27 


23.73 +0.16 


sz 


VLT/ISAAC 


2004 01/02 


> 24.3 - 24.4 


SZ 


VLT/ISAAC 


2003 04 12 


23.61 +0.15 


J 


VLT/ISAAC 



All observations of Abell 1835 should have had the same 
pointing as the archival data, namely 1' off the cluster center to 
avoid a very bright star falling into the FOV. However, a coor- 
dinate error meant that the observations were obtained at two 
pointings: the intended one and the one without the offset. The 
parts of the images that overlapped with the archival data were 
analyzed using the reference image from the archival data. The 
detection limits are given as "area 1" in Table 1. Those parts of 
the images that did not overlap with the archival data were an- 
alyzed using our observation on 18 July 2003 as the reference. 
This is "area 2" in Table 1. Because of the different photomet- 
ric depth of the reference images, we obtain slightly shallower 
detection limits for area 2. 

5. Results 

The search yielded the detection of two transient objects, one 
in each of the fields of Abell 1689 and AC114. Both objects 
were detected in the archival images and not in our new survey 
images. 

5.1. The AC 11 4 transient 

Stamps from the / and ACS F850LP-band (effectively corre- 
sponding to Gunn z-band) images of the transient in AC1 14 are 
shown in Fig. 2. The transient is seen in all archival J, H, and 
K images as a compact unresolved source. It is also present - 
barely resolved from the nearby galaxy - in the ground B and 
V band images, as well as on the ACS z-band images taken 
14 November 2004. The coordinate matching of the /-band to 
the ACS image shows that the variable object coincides exactly 
with a compact barely resolved galaxy. It is also within the po- 
sition errors of the X-ray source AC 114-2 detected by Martini 
et al. (2006) in Chandra images. The spectrum of Martini et al. 
(2006) is relatively blue, which together with the detected X-ray 
emission suggest that the object is a variable AGN. Martini et al. 
(2006) measure a redshift of ~ 0.32, which implies that the ob- 
ject is in AC 114. We also subtracted the 2 epochs of archival 
ISAAC //-band imaging obtained on 20 August 2002 (coeval 
with the /-band observations) and 17 May 2004. The subtrac- 
tion does not show a significant flux change. AGNs may vary 
slowly on time scales of years so the //-band results do not nec- 
essarily contradict the AGN hypothesis. Given all the available 
evidence, a variable AGN seems the most likely explanation for 
the transient in AC1 14. 

5.2. The Abell 1689 transient 

The transient in Abell 1689 was detected as a negative source in 
the S Z-band subtractions that used a stack of all archival data as 
a reference image. The archival SZ data were obtained in 4 runs 



between February 2003 and February 2004 (data under "point- 
ing 1" in Table 1), which made it possible to build a light curve. 
The magnitude of the transient was SZ — 23.24 + 0.08 on 9 
February 2003. It was half a magnitude fainter on 27 April 2003 
and was not detected in the combined January/February 2004 
image. The 90% detection limit of the latter image is ~ 24.15 
mag. However, the images with added artificial stars were ex- 
amined and it was noticed that if the object location is known, 
objects down to 24.3-24.4 mag could be detected visually. 

Abell 1689 has been observed on many occasions at the VLT, 
including times around the epoch at which the transient was de- 
tected. The cluster has also been observed with WFPC2 and 
ACS on HS T. All these archival data made it possible to ob- 
tain four additional measurements of the magnitude of the tran- 
sient in three bands - Bessell /, Gunn z, and the NIR / band. 
The instrumental magnitudes of the transient were measured on 
reference-subtracted images with the PSF fitting technique us- 
ing DAOPHOT in IRAF. The uncertainties were estimated by 
adding artificial stars with the same magnitude as the transient 
and measuring their magnitudes in exactly the same way as the 
transient. This procedure provides a more reliable error estimate 
than the formal output of DAOPHOT because the noise in the 
reference-subtracted images is correlated. As a by-product, we 
also obtain the accuracy with which the positions of artificial 
stars could be recovered, which we take as a measurement of the 
accuracy of the measured position of the transient. The magni- 
tudes of the transient are listed in Table 2, and Table 3 lists de- 
tails of the observations used. Stamps of the images where the 
transient was present, the references, and the difference between 
them are shown in Fig. 3. Also shown are the non-detection in 
the S Z-band observations during January /February 2004 and the 
stamp from the ACS z-band image. 

The FORS2 /-band observations were obtained in February, 
June, and July 2003 and in 2001. The deep April 2001 images 
served as a reference. The magnitude of the transient could only 
be measured in the February 2003 data and the rest were too 
shallow for the transient to be detected. The FORS2 Gunn z ob- 
servations were taken on 9, 10, 26, and 27 February 2003, but 
to increase the S/N, the images for 9, 10 and 26, 27 were com- 
bined. The ACS z-band image (obtained on 16 June 2002) was 
used as a reference and yielded very clean subtractions. For both 
epochs we measured the transient at about the same brightness of 
z — 23.9 mag. The constant brightness of the transient could be 
further verified by subtracting the two FORS2 z-images, which 
left no noticeable residuals. 

The photometric calibration of the FORS2 /-band data was 
obtained with Landolt (1992) standards observed during the 
same night. The night was photometric, but as a consistency 
check, we analyzed large number of additional FORS 1/2 /-band 
observations of Abell 1689 obtained during 2001 and 2003, most 
of which were acquisition images for spectroscopic observa- 
tions. Most of the nights were photometric, and Landolt (1992) 
standards were also observed. We used them to calibrate a num- 
ber of faint stars in the field assuming 0.04 mag per airmass 
atmospheric extinction in /. The typical RMS of the calibrated 
stars was ~ 0.02 mag, confirming that the nights, including 6 
February 2003, were photometric. 

The calibration of the FORS2 Gunn z observations was ob- 
tained through the primary white dwarf standard star GD 153 
observed on the photometric night of 10 February 2003, im- 
mediately after the observations of Abell 1689 6 . We used the 



6 Stars with calibrated z magnitudes from Hamuy et al. (2001) were 
also observed but were all saturated. 
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Fig. 3. The subtraction of the images for Abell 1689. The ACS z-band image is also shown. 



Table 3. Additional archival observations used to measure the bright- 
ness of the Abell 1689 transient, with HST/ACS data and those used as 
reference images not listed. 



Date 


Exposure Seeing 


Note 




[min] [arcsec] 




VLT/ISAAC 5 Z-band 


2003 02 09" 


159 0.52 


transient 


2003 04 27 


43 0.43 


transient 


2004 01/02 


86 0.55 


non-detection 


VLT/ISAAC 7-band 


2003 04 12 


132 0.47 


transient 


VLT/FORS2 /-band 


2003 02 06 


36 0.66 


transient 


VLT/FORS2 z-band 


2003 02 09/10 


256 0.60 


transient 


2003 02 26/27 


324 0.70 


transient 



- average of observations obtained on 5,11 and 15 February. 



synthetic VRI magnitudes of GD 153 from Holberg & Bergeron 
(2006), together with the relation between the Landolt VRI and 
the Gunn z magnitudes derived by Krisciunas et al. (2004), to 
compute the z band magnitude of GD 153, finding z = 13.77 + 
0.04. Most of the uncertainty comes from the RMS around the 
relation of Krisciunas et al. (2004). We then computed the z 
band magnitudes of the stars that were used to cross-check the / 
band photometry. We assumed 0.04 magnitudes of extinction per 
airmass for the z-band. These stars were then used to calibrate 
the photometry of the transient in the combined 9-10 and 26-27 
February 2003 images. 

The single /-band observation was obtained on 12 April 
2004. As a reference, we used a /-band image from 
VLT/HAWK-I that was obtained as part of our new program to 
search for SNe behind Abell 1689 (a sequel to the pilot ISAAC 
search we present here). ISAAC //-band data from April and 
May 2003 were also available, but were too shallow to detect the 
transient. 

The coordinates of the transient were inferred from astro- 
metric solution using coordinates of the field objects from the 
Sloan Digital Sky Survey Data Release 6. The average of the 
six measured positions is a^ooo = 13 h ll m 28*170 (+0.002) and 
£2000 = -01°19'25'.'128 (+0.012), where the uncertainties are 
the standard deviation of the mean. 



6. Discussion 

The search yielded two transient objects, one in the field of 
AO 14 and one in Abell 1689. The transient in AO 14 is very 
likely a variable AGN at the cluster redshift. Below we argue that 
the transient in Abell 1689 was probably a Type IIP SN at red- 
shift z - 0.6 and that other alternatives, including a supernova at 
the cluster redshift, are unlikely. 

Figure 4 shows a zoom of the deep ACS z-band image 
(16600 s total exposure time; the other coeval deep ACS F475W, 
F625W and F775W images show the same galaxy morphology) 
centered on the region near the position of the transient behind 
Abell 1689. Four artificial stars with magnitudes 24.32 (the tran- 
sient magnitude on the ground based z images), 25, 26, and 27 
mag were also added to demonstrate the image depth. The uncer- 
tainty of the position of the transient, the standard deviation of 
the mean, was estimated from the coordinate transformations be- 
tween the images where the transient was detected and the ACS 
image. This uncertainty is ~ 0.5 ACS pixels (1 ACS pixel = 
0.05"), which is consistent with the uncertainties of the individ- 
ual position measurements ~ 1 — 1.5 ACS pixels as computed 
by adding in quadrature the uncertainty of the positions of the 
transient as estimated from the recovered artificial stars, and the 
uncertainty of the coordinate transformations to the ACS image. 

The transient was not detected in the deep June 2002 ACS 
images to a limiting magnitude z — 26.5, but a long blazar out- 
bursts with an amplitude > 2.5 mag could, in principle, offer 
an alternative explanation of the transient. The synthetic I - SZ, 
z-SZ, and SZ-J color indices of the QSO template of Glikman 
et al. (2006) match those of the transient, if the template is 
moderately reddened with E(B - V) - 0.2 and redshifted to 
z ~ 0.6 - 0.7. Some blazars do show repeated large-amplitude 
outbursts on a time scale of months to a few years, but most 
show variability with amplitude < 1 mag. However, the archival 
HS T/WFPC2 images taken in June 1996, June 1997, and April 
2007 (the lattest coeval with our new ISAAC data), the deep 
2001 VLT, and 2002 HST/ACS observations do not show any 
object at the position of the transient. Furthermore, the closest 
galaxy has an irregular shape, but the position of the transient 
is ~ 0.5"away from its center and does not coincide with any 
of the outer spots. Along with the non-detection of any X-ray 
source close to the transient position (Martini et al., 2006), the 
variable AGN hypothesis seems quite unlilkely. 

A galactic origin for the transient is also unlikely. Any plau- 
sible candidate would be either bright enough to be detected in 
the deep ACS images (e.g. novae) or would have faded away 
faster than the transient (e.g. flaring late type stars). For example, 
WZ Sge-type dwarf novae are faint in quiescence, My - +11 
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Fig. 4. ACS z-band image of the putative host galaxy. The cross marks 
the average position of the transient from the six detections and the 
size along the x and y axes is the standard deviation of the mean. Four 
artificial stars with the indicated magnitudes were added to the image 
for comparison. 



mag (Harrison et al., 2004), and may not be detected in the deep 
ACS image if farther away than 12-13 kpc. WZ Sge stars show 
large amplitude outbursts (up to ~ 7 mag), which may bring 
them above the detection limit. However, the duration of the out- 
bursts is typically less that 80 days and the brightness decreases 
faster than observed in the Abell 1689 transient. Additional ar- 
guments against this scenario are the observed red colors of the 
transient and its / magnitude. Dwarf novae at outburst have blue 
colors in contrast to the observed red colors of the transient. 
The absolute / magnitude of WZ Sge at outburst maximum is 
~ +5 (Harrison et al., 2004; Howell et al., 2004) and the ob- 
served / =* 24. 1 would imply a distance of about 70 kpc, which 
is unlikely. 

The slow decline of the SZ light curve suggests a Type IIP 
supernova. The interaction of the expanding ejecta with circum- 
stellar material in Type Iln or possibly in other SNe types (e.g. 
SNe 2002ic and 2005gj, Hamuy et al., 2003; Aldering et al., 



2006) may also lead to slow luminosity decline. A Type la SN 
at a few months after maximum should also be considered as 
a possible candidate for the observed transient possibility. The 
brightness of prompt optical emission of gamma ray bursts typ- 
ically declines much faster than the transient in Abell 1689 and 
the most plausible conclusion is that the transient was a super- 
nova. 

Figure 5 displays a 20x20"region centered on the transient. 
At the cluster redshift, this corresponds to 60x60 kpc. None of 
the galaxies close to the transient have spectroscopic redshifts. 
The photometric redshift of the closest galaxy based on B-to-K- 
band photometry (for details see Paper II) is z = 0.60 + 0.15. 
The best-fitting template is based on single starburst evolved 
by 50 Myr and only slightly reddened. The other possible host 
galaxies, the faint objects right of the transient, and the two 
bright ellipticals have photometric redshifts between 0.4 and 0.6. 
Therefore core-collapse SN in a cluster or foreground galaxy is 
unlikely. Intercluster, hostless SN la at about 80-100 days past 
maximum light could match the observed brightness of the tran- 
sient. However, the late-time (> 150 days past maximum) bright- 
ness of SNe la typically declines by ~ 1.5 mag per 100 days and 
even faster between 70 and 150 days. Thus the transient should 
have declined by about 1 mag over the 77 days between the two 
SZ observations, which is twice more than observed. Thus, we 
find it quite unlikely that the transient could be a supernova at 
the cluster redshift. 

To further investigate the nature of the transient we fit- 
ted its measured magnitudes with several templates based on 
the well-observed Type IIP SNe 2004et (Sahu et al., 2006), 
1999em (Hamuy et al., 2001; Krisciunas et al., 2009), 2003hn 
(Krisciunas et al., 2009), and 2001cy (Poznanski et al., 2009). 
The SNe were carefully selected to sample all light curve mor- 
phologies of Type IIP SNe. We also used 2005gj (Prieto et al., 

2007) , which had a slow decline in brightness after maximum 
light because of interaction with the circumstellar medium. The 
three SNe with NIR photometry (2003hn, 1999em, and 2005gj) 
also allowed us to test redshifts below 0.4. Because only 1999em 
have (few early) NIR spectra, the NIR magnitudes of these SNe 
were fitted with a cubic spline and the resulting fits were used to 
compute the cross-filter /f-corrections (Kim et al., 1996). 

The templates were used to perform lightcurve fits on a grid 
of redshifts and reddening parameters. Since we only have six 
measurements, we tried three fixed values of Ry and only al- 
lowed E(B - V) to vary, along with 3 other lightcurve parame- 
ters: redshift, time of maximum, and peak magnitude. The best 
fit, shown in Fig. 6, was obtained with the template based on 
SN 200 Icy redshifted to z = 0.59 + 0.05 and reddened by non- 
standard dust in the host galaxy with Eb-v = 0.85 and Ry = 1.5. 
It should be noted that the uncertainty of the redshift estimated 
from the light curve fit is the formal error from the fit and may be 
somewhat underestimated. The lower panel of Fig. 6 show the^ 2 
as a function of the redshift for the three best-fitting templates, 
all for R v = 1.5, as well as SN 2001cy for the standard Milky 
Way Ry = 3.1. It is interesting to note that the best fits were ob- 
tained for low Ry — 1.5 and redshift around 0.6 and that lower 
redshifts result in much poorer fits. SNe 1999em and 2004et do 
not fit the data well because their R and / brightness are either 
constant or increases on the plateau phases, while the transient 
SZ magnitude decreases. 

The low best-fit value of Ry deserves special comment. 
Recent studies of individual nearby SNe la and statistical anal- 
ysis of large samples of SNe la (e.g., Elias-Rosa et al., 2006, 
2008; Astier et al., 2006; Nobili & Goobar, 2008) and now also 
for Type IIP SNe (Poznanski et al., 2009; Olivares, 2008) find 
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SN2001cy: % =3.5 E(B-V)=0.85 R v =1.5 



Fig. 5. Zoom of the ACS z-band image around the transient in Abell 
1689. The image size is 20x20" or 60x60 kpc at the cluster redshift. 
The position of the transient is marked and the estimated photometric 
redshift of some galaxies is also shown. 



that dust in SN host galaxies has Ry significantly lower than the 
average Milky Way value of ~ 3.1. It is possible that the dust 
that obscured the transient had properties corresponding to low 
Ry, possibly as a result of multiple scattering by circumstellar 
dust clouds (Goobar, 2008). It should be noted that in the fitting 
process, we used the Milky Way extinction law parametrization 
by Fitzpatrick (1999), which was derived using observations of 
stars obscured by dust with Ry between 2.3 and 5.3 7 . Any use of 
this parametrization to compute the extinction laws for Ry < 2.3 
thus involves extrapolation, which may be inaccurate, especially 
if Ry <k 2.3. A more reliable approach may be to use a physical 
model of the interstellar dust (Weingartner & Draine, 2001), but 
this is beyond the scope of this paper. 

To estimate the absolute magnitude of the transient, we cor- 
rect the peak magnitude that we estimated from the light curve 
fit for the dust extinction and gravitational magnification of the 
cluster. For this we used the accurate parametric mass distri- 
bution model of Abell 1689 from Limousin et al. (2007). The 
model is based on strong lensing analysis and uses 34 multi- 
ply imaged systems, 24 of them with newly determined spec- 
troscopic redshifts. As discussed by Limousin et al. (2007), the 
large number of constraints used to derive the mass model make 
Abell 1689 the most reliably reconstructed cluster to date. For 
more details on the mass model of Abell 1689 see Limousin et al. 
(2007). 

The cluster magnification at the position of the transient for 
redshift z = 0.6 is 1.36 +0.03 mag, where the uncertainty is es- 
timated from the uncertainties of the mass model. For redshift 
z = 0.45 and z - 0.75, the magnifications are 1.13+0.02 and 
1 .60+0.04, respectively. The uncertainty from the mass model is 
much smaller than the uncertainty stemming from the photomet- 
ric redshift estimation and we adopt a magnification of 1.4+0.3 
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Fig. 6. Upper: Best fit of the transient magnitudes with a light curve 



template based on SN 2001cy. Lower: x vs - 
collapse SN light curve templates. 
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7 The other widely used parametrization of Cardelli et al. (1989) was 
obtained with stars obscured by dust with R v between 2.75 and 5.3 



mag to estimate the absolute magnitude of the transient. Taking 
the dust extinction into account, we find absolute V magnitude 
of the transient My ~ -17.6 + 0.3, which is normal for a Type 
IIP SN (e.g., Richardson et al., 2002). The transient redshift in- 
ferred from the template fit matches the photometric redshift of 
the nearest galaxy very well, which is the most likely host of the 
SN. With its young stellar population inferred from the best pho- 
tometric redshift fitting galaxy template and the irregular mor- 
phology, this galaxy is likely to produce core-collapse SNe. 



7. Conclusions 

We have presented the first dedicated rolling search for gravita- 
tionally magnified SNe behind galaxy clusters. The search was 
conducted in the near infrared wavelengths using the ISAAC in- 
strument at the VLT. The massive galaxy clusters Abell 1689, 
Abell 1835 and AC 114 were observed for a total of 20 hours, 
split into individual epochs of 2 hours exposure time each sepa- 
rated by approximately one month. Image subtraction was used 
to search for transient objects with light curves consistent with 
supernovae. The analysis of the combined data set of new obser- 
vations, plus 20 more hours of archival data, yielded the detec- 
tion of two transient objects, both found in archival images of 
Abell 1689 and AC114. The transient in AC114 coincides with 
an X-ray source and is probably a variable AGN at the clus- 
ter redshift. The transient in Abell 1689 was found with an SZ 
magnitude of 23.24 and ~ 0.5"away from a galaxy with photo- 
metric redshift z ga i = 0.6 + 0.15. The light curves and the col- 
ors of the transient are consistent with it being a reddened Type 
IIP supernova at redshift z = 0.59 + 0.10. The lensing model 
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of Abell 1689 predicts 1.4±0.3 mag of magnification at the po- 
sition of the transient, making it the most magnified supernova 
ever found and only the second supernova found behind a galaxy 
cluster. Without the magnification, this SN would not have been 
detected in our NIR search. 

The limiting magnitude of the individual epochs was esti- 
mated by adding artificial stars to the subtracted images. Most of 
the epochs reach 90% detection efficiency at S Z( J) =* 23.8-24.0 
mag (Vega). The detailed analysis presented in the accompany- 
ing Paper II indicates that our survey should have detected be- 
tween 0.8 and 1.6 SNe with the peak probability corresponding 
to Type IIP SNe at a redshift of z ~ 0.7. It is encouraging that 
this is what we actually found. 

Our pilot survey has demonstrated the feasibility to find dis- 
tant SNe gravitationally magnified by massive galaxy clusters. 
The survey reached a detection threshold of ~ 23.8 - 24.0 mag. 
A deeper survey of several clusters with an NIR camera with a 
larger field-of-view should yield an order of magnitude increase 
in the number of detected SNe, a significant fraction of which 
will be at redshifts z > 1.5 (Paper II). 
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